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ABSTRACT 

In  this  paper,  two-dimensional  equation,  subject  to  Dirichlet  initial  boundary  conditions  is  presented,  and  a 
Spline  Collocation  Method  is  utilized  for  solving  the  problem.  Also,  Spline  provides  continuous  solution,  in  contrast  to 
another  method,  which  only  provides  discrete  approximations.  It  is  found  that  this  method  is  a powerful  mathematical 
tool  and  can  be  applied  to  a large  class  of  linear  and  nonlinear  problems  in  different  fields  of  science  and  technology. 
Numerical  results  obtained  by  the  present  method  are  in  good  agreement  with  the  analytical  solutions  available  in  the 
literature. 
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INTRODUCTION 

Recently,  new  numerical  methods  have  gained  the  interest  of  researchers  for  finding  numerical  solutions 
to  partial  differential  equations.  This  interest  was  driven  by  the  need  for  applications,  both  in  industry  and 
sciences.  Theory  and  numerical  methods  for  solving  initial  boundary  value  problems  were  investigated  by  many 
researchers.  In  the  last  decade,  there  has  been  a growing  interest  in  the  new  numerical  techniques  for  linear  and 
non-linear  initial  boundary  value  problems.  In  the  field  of  partial  differential  equations,  boundary  value  problem  is 
a differential  equation,  together  with  a set  of  additional  restraints,  called  the  boundary  conditions.  A solution  to  a 
boundary  value  problem  is  a solution  to  the  differential  equation,  which  also  satisfies  the  boundary  conditions. 
Boundary  value  problems  arise  in  engineering,  applied  mathematics  and  several  branches. 

However,  it  is  usually  difficult  to  obtain  closed-form  solutions  for  boundary  value  problems,  especially 
for  non-linear  boundary  value  problems.  In  most  cases,  only  approximate  solutions  (either  numerical  solution  or 
analytical  solutions)  can  be  expected.  Some  numerical  methods  such  as  finite  difference  method,  finite  element 
method,  spline  approximation  method,  Jacobi’s  Method  and  SOR  method  have  been  developed  for  obtaining 
numerical  solutions  to  boundary  value  problems. 

In  this  paper,  we  apply  a spline  collocation  method  approach,  to  obtain  a solution  to  the  wider  class  of 
non-linear  systems  of  boundary  value  problems.  The  spline  method  approach  is  widely  utilized  for  the  numerical 
solution  of  non-linear  problems  arising  in  real  world  applications. 

The  governing  equations  here  are  highly  partial  differential  equations,  which  are  solved  by  using  the 
spline  collocation  method.  In  this  way,  the  paper  has  been  organized  as  follows.  In  section  2,  we  use  the  spline 
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collocation  method,  in  section  3 spline  solution  of  explicit  method  and  implicit  method.  The  conclusions  are  summarized  in 
section  4. 


SPLINE  COLLOCATION  METHOD 

The  general  form  of  hyperbolic  PDE  with  one  space  variables  x 
and  time  variable  t is  given  by 


with  a Dirichilet  boundary  conditions,  namely 

u(0,  y,  t)  = u(x,  0,t)  = 

and  two  initial  conditions  at  t = 0 (Cauchy  conditions) 


(i) 


(ii) 


u(x,y,0)  = f(x,y) 

ut  (x,  y,  0)  = g(x,  y).  (iii) 

In  equation  (i),  c2  is  a constant  term;  it  depends  upon  some  physical  quantities  in  case  of  different  types  of 
problems. 

Divide  the  region  0<  x < a into  say  n sun-intervals  each  of  width  Ax  such  that  nAx  = a,  0<  y < b into  say  m sun- 
intervals  each  of  width  Ay  such  that  nAy  = b. 

The  subscript  k denotes  the  time  and  (I,  j)  denote  the  position.  For  simplicity,  consider  square  region,  i.e.  a = b = 
L the  length  of  the  region, 

m = n and  Ax  = Ay  = h. 

The  points  of  subdivisions  are  (xi?  yi),  i,j  = 0(1  )n. 

For  explicit  scheme,  the  formula  is  obtained  in  a similar  manner  as  discussed  in  section  (3.9),  by  replacing  the 
right  side  by  twice,  the  second  derivative  of  cubic  spline  S (x)  at  the  (i,  j) th  mesh  point  I.  e.  2S"  (xi;j). 

Let  ui,j  we  get 

(Ui-lj,  k+l+  4uy9k+i+  ui+i9j9k+i) 

= (2+12r2)  ui+Uk+(8-24r2)  uy,k+(2+12r2)  u,.,j,k 

-[  Ui+ij,k.i+4uyjk_i+  Ui_ij9k_i]  (iv) 

Where  r = cAt/h  i,j  = l(l)n-l 

The  values  of  u^.i,  i,j  = 1(1)  n -1  are  obtained  from  following  relation 

Uij,-i  = Uij.,1  - 2g(Xi,yj)At  (v) 

The  equation  (iv)  is  known  as  a cubic  spline  explicit  formula  to  solve  hyperbolic  type  PDE  of  the  form  given  by 


Impact  Factor  ( JCC ):  6.8765 


NAAS  Rating:  3.11 


Vibrating  Membrane  Problem  Solved  using  Spline  Collocation  Method  with  Dirichlet  Conditions 


149 


the  equation  (i). 

The  values  obtained  from  equation  (v)  using  initial  and  boundary  conditions,  the  equation  (iv)  gives  a system  of 
(n-1)  x (n-1)  simultaneous  equations  in  (n-1)  x (n-1)  unknowns. 

After  calculating  the  results  for  k = 0;  the  results  for  k = 1,2,....  are  obtained  in  a similar  manner. 

r < (1/2)  is  the  required  condition  for  convergence  and  stability  of  this  explicit  method. 

Likewise,  implicit  scheme  for  obtaining  the  solution  of  parabolic  PDE  with  one  space  variable,  we  have  an 
implicit  scheme  to  solve  hyperbolic  differential  equation  with  two  space  variables. 

This  implicit  scheme  is  unconditionally  stable. 

The  cubic  spline  implicit  formula  to  solve  equation  (i)  is  as  follows,  by  replacing  the  right  side  of  equation  (i) 

(l-6r2)ui+1J>k+1+  (4+12r2)uiJ>k+i+  (l-6r2)Uj_1>j>k+1 

= (6r2-l)  ui+iJ>k.1+(4+12r2)  uiJ>k.|+(6r2-l)  Ui.1J>k.i 

+2[  ui+i J>k+4ui!j>k+  uMJ,k]  (vi) 

Where  r = cAt/h  i,j  = l(l)n-l 

Due  to  similar  discussions  as  above,  also  equation  (vi)  gives  (n-1)  x (n-1)  simultaneous  equations  in  (n-1)  x (n-1) 
are  unknowns. 

After  calculating  results  for  (k+l)th  time  level,  the  results  for  (k+2)th  time  level  are  obtained  in  similar  manner  by 
repeating  the  process. 

The  problem  of  vibrating  membrane  is  well-known  to  researchers.  Assume  that  membrane  is  tightly  stretched  and 
homogeneous  i.e.  its  mass  per  unit  area  is  constant,  it  is  perfectly  flexible  and  is  so  thin  that,  it  offers  no  resistance  for 
bending. 

The  membrane  is  stretched  and  then  fixed  along  its  entire  boundary  in  the  xy  - plane,  and  the  tension  T caused  by 
stretching  the  membrane  is  the  same  at  every  point  in  all  directions  and  does  not  change  during  the  motion. 

The  deflections  u (x,  y,  t)  of  the  membrane  during  the  motion  are  small  compared  to  the  size  of  the  membranes. 

Consider  the  forces  acting  on  a small  portion  of  the  membrane,  as  shown  in  below  figure. 

Since  the  deflection  of  the  membrane  and  the  angles  of  inclination  are  small,  the  sides  of  the  portion  may  be  taken 
approximately  equal  to  Ax  and  Ay. 

The  tension  T is  the  force  per  unit  length.  Forces  acting  on  the  edges  are  TAx  and  TAy,  which  could  be  taken 
tangent  to  the  membrane,  since  the  membrane  is  perfectly  elastic. 

Let  the  force  TAy  make  angles  A and  B with  the  horizontal  on  the  opposite  edges  of  the  membrane,  since  the 
membrane  is  perfectly  elastic. 

The  resultant  vertical  component  of  force  due  to  TAy  is  therefore, 

(TAy)  sinB  - (TAy)  sin  A 
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(TAy)tanB-(TAy)tanA 
= TAy{  (5u/5x)x+Ax-(5u/5x)x) 

= TAy{(S2u/Sx2}Ax 

Up  to  a first  order  approximation,  note  that  sine  have  been  replaced  by  tangent,  because,  the  angles  A and  B are 

small. 

Similarly,  the  forces  TAx  acting  on  the  edges  of  length  Ax  can  be  shown  to  have  the  vertical  component. 
TAxKdVdy^Ay 

If  m is  the  mass  per  unit  area  of  the  membrane,  by  Newton’s  second  law  of  motion 
(mAxAy)  d2u/dt2  = T {(d2u/dx2+d2u/dy2}AxAy 

d2u/dt2  = c2{(d2u/dx2+d2u/dy2}  where  c2  = T/m  .(vii) 

The  above  equation  (vii)  is  hyperbolic  PDE  in  two  space  variables  x and  y,  and  time  variable  t. 

The  solution  for  the  case  of  a rectangular  membrane  is  given  below  by  cubic  spline  explicit,  as  well  as  an  implicit 
method,  using  the  following  initial  and  boundary  conditions. 

Dirichilet  boundary  conditions: 

u(0,y,t)  = 0 0 < y < t 

u(x,0,t)  = 0 0 < x < t (viii) 

Initial  conditions  at  t = 0 (Cauchy  condition) 

u(x,y,0)  = f(x,y)  = simix.simry  ; 0 < x,y  < 1 (ix) 

u(x,y,0)  = 0 

SPLINE  SOLUTIONS  WITH  EXPLICIT  METHOD 

Using  initial  and  boundary  conditions  are  described  by  equations  (viii)  and  (ix),  the  solution  of  equation  (vii)  is 
obtained  through  explicit  formula  given  by  the  equation  (iv)  as  follows. 

Let  the  region  0 < x < 1 and  0 <y<l  be  divided  into  10  sub -intervals,  each  of  length  h = 0.1.  Let  c = 0.1  and  At  = 
0.01,  these  gives  r = 0.01.  After  calculating  the  values  of  uis  j_i,  i,  j = 1(1)  n form  equation  (v),  using  initial  condition, 
substitute  the  values  of  ui;  j_i,  2 + 12r2  and  8 - 24  r2  into  the  equation  (iv)  one  gets, 

Lor  k = 0 

Lor  j = 1 

i = 1 Uq,  i;  i + 4ui  11  + U2;  i,  i = 0.563599 

Since  u0,  i,  i = 0 

4u}> ! i + u2>  i,  i = 0.563599 

i = 2 ui,  11  + 4u2,  i,  i + U3  i i=  1.072029 


Impact  Factor  (JCC):  6.8765 


NAAS  Rating:  3.11 


Vibrating  Membrane  Problem  Solved  using  Spline  Collocation  Method  with  Dirichlet  Conditions 


151 


i = 3 U2;  i5 1 + 4u3  ii  + U4  i i=  1.475521 
i = 4 U3  i,  1 + 4u4  11  + U5  ii=  1.734579 
i = 5 U4  1 1 + 4u5,  1, 1 + U6  11=  1.873844 
Since  14, 1, 1 = U6,  it  1 (Due  to  symmetry  of  the  problem) 

2u4?  if ! + 4u5?  ] 1 = 1.873844 

In  a similar  manner,  calculate  the  equation  j = 1 (1)9  one  gets  9x9  simultaneous  linear  equations  in  9 x 9 
unknowns. 

This  system  of  equations  can  be  solved  by  any  well-known  method.  Once  the  values  of  u are  known  at  the  first 
level  of  time,  the  process  can  be  repeated  for  second  time  level  and  so  on.  The  results  obtained  by  explicit  method  are 
given  as  follows.  Due  to  the  symmetry  of  the  problem,  results  are  given  for  0 < y < 0.5,  0 < x < 0.5  in  the  table  (a)  and  (b). 

The  results  obtained  at  y = 0.1  are  plotted  in  the  figures  (P). 


Table  1:  Velocity  Distribution  in  Rectangular  Vibrating  Membrane  through 
Cubic  Spline  Explicit  Method  u— ► at  t=0.01 


x\y 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0 

0 

0 

0 

0 

0 

0 

0.1 

0 

0.095491 

0.181634 

0.249998 

0.293890 

0.309014 

0.2 

0 

0.181634 

0.345488 

0.475523 

0.559012 

0.587779 

0.3 

0 

0.249998 

0.475523 

0.654502 

0.769413 

0.809009 

0.4 

0 

0.293890 

0.559012 

0.769413 

0.904499 

0.951047 

0.5 

0 

0.309014 

0.587779 

0.809009 

0.951047 

0.999990 

Table  2:  Velocity  Distribution  in  Rectangular  Vibrating  Membrane  through 
Cubic  Spline  Explicit  method  u — > at  t=0.03 


x\y 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0 

0 

0 

0 

0 

0 

0 

0.1 

0 

0.095483 

0.181619 

0.249978 

0.293867 

0.308989 

0.2 

0 

0.181619 

0.345461 

0.475485 

0.558967 

0.587732 

0.3 

0 

0.249978 

0.475485 

0.654450 

0.769352 

0.808944 

0.4 

0 

0.293867 

0.558967 

0.769352 

0.904427 

0.950917 

0.5 

0 

0.308989 

0.587732 

0.808944 

0.950971 

0.999910 

Figure  1 
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SPLINE  SOLUTIONS  WITH  IMPLICIT  METHOD 


In  similar  manner  discussed  as  earlier,  using  initial  and  boundary  condition  at  At  = 0.01,  r = 0.01,  the  solution  of 
equation  (i)  is  obtained  by  using  cubic  spline  implicit  formula  given  by  equation  (vi)  as  follows. 

For  k = 0 


For  j = 1 

i = 1 (4.0012)  ult  i,  i + (0.9994)  u2,  i,  i = 0.563602 

I = 2 (0.9994)  ulf  i,  i + (4.0012)  u2f  i,  i + (0.9994)  u3,  i,  i = 1.072034 

i = 3 (0.9994)  u2  lf  x + (4.0012)  u3,  i,  i + (0.9994)  u4,  i,  i = 1.475528 

I = 4 (0.9994)  u3  lf  i + (4.0012)  u4,  i,  i + (0.9994)  u5>  lt  r = 1.734587 

i = 5 (1.9988)  u4  !,  i + (4.0012)  u5,  i,  i = 1.823853 

One  gets  system  of  9x9  simultaneous  linear  equations  in  9x9  unknowns. 

Once  the  values  of  u are  obtained  at  the  first  time  level,  the  values  for  second  time  level  are  obtained  in  similar 
manner  by  repeating  the  process. 

Results  obtained  by  cubic  spline  implicit  method  are  given  at  different  times  in  the  below  tables. 

The  results  are  plotted  in  the  below  figures  at  y = 0. 1 . 


Table  3:  Velocity  Distribution  in  Rectangular  Vibrating  Membrane  through 
Cubic  Spline  Implicit  Method  u — > at  t = 0.01 


x\y 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0 

0 

0 

0 

0 

0 

0 

0.1 

0 

0.095491 

0.181634 

0.249998 

0.293890 

0.309014 

0.2 

0 

0.181634 

0.345488 

0.475523 

0.559012 

0.587779 

0.3 

0 

0.249998 

0.475523 

0.654502 

0.769413 

0.809009 

0.4 

0 

0.293890 

0.559012 

0.769413 

0.904499 

0.951047 

0.5 

0 

0.309014 

0.587779 

0.809009 

0.951047 

0.999990 

Table  4:  Velocity  Distribution  in  Rectangular  Vibrating  Membrane 
through  Cubic  Spline  Implicit  Method  u — >at  t = 0.03 


x\y 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0 

0 

0 

0 

0 

0 

0 

0.1 

0 

0.095491 

0.181634 

0.249998 

0.293890 

0.309014 

0.2 

0 

0.181634 

0.345488 

0.475523 

0.559012 

0.587779 

0.3 

0 

0.249998 

0.475523 

0.654502 

0.769413 

0.809009 

0.4 

0 

0.293890 

0.559012 

0.769413 

0.904499 

0.951047 

0.5 

0 

0.309014 

0.587779 

0.809009 

0.951047 

0.999990 
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Figure  2 


DISCUSSIONS  OF  RESULTS 


The  velocity  distribution  obtained  from  the  solution  of  equation  (i)  by  cubic  spline  explicit  as  well  as  implicit 
methods  are  compared  with  exact  solutions  at  y = 0.1  for  t = 0.03  in  the  above  table. 

From  the  table,  it  is  clear  that  the  solutions  obtained  by  both  these  methods  are  accurate,  up  to  five  digits  of 
decimal  places. 

The  figure  indicates  analysis,  which  compares  the  exact  solutions  with  spline  solutions  obtained  by  both  the 
methods,  explicit  as  well  as  implicit  at  t = 0.03,  clearly  revealing  that  the  spline  solutions  are  quite  accurate  and  reliable. 
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